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ABSTRACT 

The dynamics of artificial asteroids on the Trojan-hke orbits around Neptune is in- 
vestigated in this paper. We describe the dependence of the orbital stability on the 
initial semimajor axis a and inclination i by constructing a dynamical map on the 
(a, i)— plane. Rich details are revealed in the dynamical map, especially a unstable 
gap at i = 45° is determined and the mechanism triggering chaos in this region is 
figured out. Our investigation can be used to guide the observations. 
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1 INTRODUCTION 

In the restricted three-body model consisting of the Sun, a 
planet and an asteroid, the equilateral triangular Lagrange 
equilibrium points (L4 and Ls) are stable for all planets in 
our Solar system. Asteroids in the vicinities of L4 and L5 
are called Trojans after the the group of asteroids found 
around Jupiter's Lagrange point. Trojan asteroids of Mars 
and Earth have also been observed while the Trojan-type 
orbits of Saturn and Uranus have been proved unstable 
due to the perturbations from other planets. As for Nep- 
tune, 6 Trojan-like asteroids were discovered in recent years 
(lAU: Minor Planet Center, http://www.cfa.harvard.edu/ 
iau/Iists/NeptuneTrojans.html). We list their orbital prop- 

^ erties in Table 1. The long-term orbital stability of these 
cisteroids has been studied and verified in different papers, 

' e.g. (Marzari et al. 2003; Brasser et al. 2004b; Li et al. 2007). 
There could be much more Trojan-type asteroids sharing the 
orbit with Neptune than with Jupiter, both in the sense of 
number and total mass (Sheppard & Trujillo 2006). There- 
fore it is worth to investigate the stable region in the whole 
parameter space. This topic can be found in several papers 
(Nesvorny & Dones 2002; Dvorak et al. 2007). 

Since the L5 point of Neptune is nowadays in the di- 
rection of the Galaxy center and not suitable for asteroid 
observing, all the asteroids in Table 1 are around the L4 
point. There are reports that the shape and size of the sta- 
ble regions around L4 and L5 points are different from each 
other (Holman & Wisdom 1993), but further analysis prove 
that this asymmetry is no more than an artificial effect of 
asymmetric initial conditions (Nesvorny & Dones 2002; Dvo- 
rak et al. 2008). Thus, it is reasonable to study only one of 
the Lagrange points and expect the other one has the same 
dynamical behavior. 

All the Trojans in Table 1 are on the near-circular orbits 
(small eccentricities) and two of them have high inclination 
values. The origin of the high-inclined orbit is an interesting 



Table 1. Orbits of Neptune's Trojans. The mean anomaly M 
is given at epoch TD=20080514. The perihelion argument u), as- 
cending node n and inclination i are in degree (J2000.0). 



Designation 


M 


iU 




i 


e 


a(AU) 


2001 QR322 


58.89 


158.5 


151.6 


1.3 


0.031 


30.262 


2004 UPlO 


339.30 


359.3 


34.8 


1.4 


0.027 


30.171 


2005 TN53 


284.82 


86.7 


9.3 


25.0 


0.064 


30.143 


2005 T074 


265.46 


304.0 


169.4 


5.3 


0.052 


30.151 


2006 RJ103 


231.55 


32.7 


120.8 


8.2 


0.027 


30.036 


2007 VL305 


351.72 


215.1 


188.6 


28.1 


0.062 


30.007 



topic (Li et al. 2007) , but in this paper, we will discuss only 
the stability of inclined orbits and try to find out the possible 
region (in dynamical sense) where the potential Neptune 
Trojans could survive for long time. 



2 DYNAMICAL MAP 

To investigate the efi'ects of inclination on the stability of 
Trojans, we numerically simulate the evolutions of thou- 
sands of test particles on the Trojan-like orbits. The dynam- 
ical system consists of the Sun, four jovian planets (Jupiter, 
Saturn, Uranus and Neptune) and the massless test parti- 
cles. For each set of initial conditions, a specific inclination 
value is given and 101 artificial Trojans are initialized around 
the L5 point of Neptune as follow: Their eccentricities eo, 
ascending nodes fio and mean anomalies Mq are exactly the 
same as the ones of Neptune, but the perihelion arguments 
differs 60° from the one of Neptune. Because the Trojans 
share the same orbit with the planet, they are in the 1:1 
mean motion resonance with the planet. The critical argu- 
ment of this resonance is ct = A — Aat where \ = ijj-\-Q.-\- M \s 
the mean longitude and the subscript 'N' denotes Neptune. 
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In this paper, we always study the case of the traihng La- 
grange point (is), so that ctq = —60°. The scrniniajor axes 
ao of test particles are from 29.9 AU to 30.5 AU with an in- 
crement 0.006 AU (the osculating semimajor axis of Neptune 
is 30.14 AU at the starting of simulation). Finally, wc vary 
their inclinations from 0° to 70° with an increment of 1.25°. 
The systems are then integrated to 2.687 x lO'^ yrs with a 
Lie- integrator (Hanslmeier & Dvorak 1984). An on-line low- 
pass digital filter is applied to filter the high-frequency terms 
in the output and reduce the final data size. Wo apply spec- 
tral analysis to the final data and use the spectral number to 
indicate the regularity of an orbit. In principle, the spectral 
number (hereafter SN) is the number of peaks in a power 
spectrum which are higher than a specific threshold. For de- 
tails of calculating SN, see for example (Ferraa-Mello et al. 
2005). 

Figure 1 shows the dynamical map on the initial (oo, io) 
plane. The grey depth indicates the SN of the resonant ar- 
gument a. The spectral number is forced to be 100 if it is 
greater than 100, and those orbits with averaged semimajor 
axis a ^ [29.9, 30.5] (AU) are given a spectral number of 110 
(dashed area in Figure 1) since they are surely not inside the 
1:1 resonance. Those orbits with small SNs are dominated by 
few dominating frequencies thus are more regular while the 
bigger SNs indicate strong noises in the motion and chaotic 
orbits. To verify the reliability of this regularity indicator, 
we also integrate hundreds of orbits to the solar system age 
(4.5G years) using the hybrid symplectic integrator in the 
Mercury 6 integrator package (Chambers 1999). The com- 
parison between the SNs derived from our 2.687 x 10^ yrs 
integrations and the results from the 4.5G years integrations 
done by the Mercury6 show that orbits with small SNs are 
generally regular and survive on the Trojan-like orbits in 
4.5G years while orbits with SNs higher than ~ 60 will es- 
cape from the resonant region and/or be ejected by (collide 
with) the planets or the Sun. 

There are some interesting features in Figure 1 deserv- 
ing a description. 

1) The center of the resonant orbits is at ~ 30.218 AU 
in terms of the initial semimajor axis while the mean semi- 
major axis is ~ 30.11 AU. The stable region is distributed 



symmetrically with respect to this center. Away from the 
center, the libration amplitude of the resonant argument Act 
increases as shown in Figure 2 and the Trojans run on the so- 
called 'tadpole' orbit. It is reported that no Neptune's Tro- 
jan can survive with Aa > 60° (Holnian & Wisdom 1993; 
Nesvorny & Dones 2002). Comparing Figure 1 and 2, we 
can derive the same conclusion. The inner and outer edges 
of the resonant region seen in Figure 1 oo ~ 30.04 AU and 
■~ 30.40 AU are defined by the overlapping of the secondary 
resonances. 

2) The stable orbits with initial inclination as high as 
io = 60° exist as Figure 1 indicates. This upper limit in fact 
can be found in a restricted three-body model in which it's 
61.7° (Brasscr et al. 2004a). Different values (35° and 70°) 
of such a threshold (Nesvorny & Dones 2002; Dvorak et al. 
2007) probably are due to the specific initial conditions of 
the orbits, the stability criteria and the inadequate sample 
Trojans in their simulations. 

3) The most distinguishable feature in Figure 1 is a 
gap locating at io ~ 45°. It separates the resonant region 
into two disconnected parts. Therefore we may expect to 
find two primordial Trojan groups in future, because this 
gap prohibits Trojans in one region from entering the other 
through secular diffusion. 

4) There are rich details in the resonant region. For 
example, two less regular regions at low inclination and with 
ao ~ 30.11 and ~ 30.31 AU are visible in Figure 1. They arise 
from the secular resonance via, that is, the nodal frequency 
of the Trojan in these regions are very close to Neptune's 
nodal frequency, fl ^ Qn- Another noteworthy structure is 
an arc of irregular motion extending from (oo = 30.04 AU, 
io — 0°) to (ao — 30.15 AU, io = 23°), and the symmetrical 
arc on the right side also. Although the mechanism behind 
is not understood very well, the arc is reflected in Figure 2 
again as a valley of large libration amplitude. 

Up to now all the initial eccentricities of test Trojans 
arc set to be the same as Neptune (eo — 0.006). In the 
dynamical evolution, the eccentricities of stable orbits are 
kept small. In fact for most of stable orbits with io < 45° , 
the maximum eccentricity is smaller than 0.05. Only some 
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Figure 3. A typical orbit in the unstable gap in Figure 1. 
The initial orbital elements of this orbit are ao = 30.218 AU, 
eo = 0.006,40 = 46.25°. Prom top to bottom, the panel shows 
the evolution of the semimajor axis, eccentricity, inclination, per- 
ihelion argument, resonant argument and the difference between 
the perihelion longitudes of the Trojan and of Neptune. The in- 
tegration of this orbit was terminated at 2.45 x lO'^ yrs when this 
Trojan was ejected by a close encounter with Uranus. 

Trojans with io > 45° may be on eccentric orbits, but the 
eccentricities are still limited by e ~ 0.1. 



3 MOTION IN THE GAP 

The most prominent trait in the dynamical map (Figure 1) 
is the gap at io ~ 45° . We applied the frequency analysis to 
orbits starting from the gap and have figured out the mech- 
anism causing the chaotic motion. In Figure 3, we illustrate 
a typical orbit in the gap and use this as an example to 
explain how an orbit in the gap evolves. 

As shown in Figure 3, the resonant argument a = X — 
\n librates with a small amplitude around the Lagrange 
point Ls with cr € (-65°, -52°) before 1.57 x 10^ yrs. This 
libration implies that the Trojan is in the 1:1 mean motion 
resonance with Neptune. Thanks to the protection of the 
resonance, the evolutions of other orbital elements during 
this period are also regular. For example the semimajor axis 
is nearly constant and the inclination variation around ~ 44° 
with an amplitude of only ~ 4°. But there is one exception, 
the eccentricity of the Trojan keeps increasing during this 
period and, the eccentricity reaches e = 0.355 at T = 1.57 x 
lO'^ yrs. We know that the secular resonance related to the 



precession of the perihelion may drive the eccentricity up 
(Murray & Dermott 1999). We check the frequency of the 
perihelion longitudes of the Trojan and the planets in the 
system, the result show that the secular resonance vs (ri? « 
lijjv) is responsible for the eccentricity increasing. A proof of 
this secular resonance is clearly shown in the bottom panel 
in Figure 3 where w — wn librates around ~ 270°. 

This high eccentricity makes the Trojan's perihelion dis- 
tance q — a{l — e) « 19.4 AU, which means the Trojan 
may cross the orbit of Uranus {au = 19.2 AU). But its high 
inclination makes the probability of close encounter with 
Uranus very small thus the orbit can be still safe. How- 
ever, the high eccentricity also makes another secular reso- 
nance possible, the Kozai resonance (Kozai 1962), in which 
the perihelion argument w librates while the eccentricity 
and inclination undergo variations such that the quantity 
Hk = Vl — cosi remains constant. These can be found 
after T = 1.57 x 10^ yrs in the 2nd, 3rd and 4th panel of 
Figure 3. 

In Kozai resonance, when the eccentricity increases the 
inclination decreases. Consequently the probability of close 
encounter with Uranus or other planets is enhanced signif- 
icantly, and such close encounters make the Trojan's orbit 
unstable, as Figure 3 shows. 



4 CONCLUSION 

The dependence of the stability of Neptunian Trojans on 
their inclinations is investigated and shown by a dynami- 
cal map. A gap of unstable orbits with initial inclination of 
~ 45° has been discovered. The mechanism responsible for 
this unstable gap, a combined effects from the secular 
resonance and the Kozai resonance has been figured out. 
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